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^— ( Abstract. This paper considers the problem of finding a quickest path between two points 

7—i in the Euclidean plane in the presence of a transportation network. A transportation network 

consists of a planar network where each road (edge) has an individual speed. A traveller may 
enter and exit the network at any point on the roads. Along any road the traveller moves with 
a fixed speed depending on the road, and outside the network the traveller moves at unit speed 
in any direction. 

We give an exact algorithm for the basic version of the problem: given a transportation network 
of total complexity n in the Euclidean plane, a source point s and a destination point t, find a 

^\ quickest path between s and t. We also show how the transportation network can be preprocessed 

in time 0(n^ log^^) into a data structure of size 0{n^) such that a (1 + £)-approximate quickest 

path cost queries between any two points in the plane can be answered in time 0(l/e* logn). 

o 

1 Introduction 

O Transportation networks are a natural part of our infrastructure. We use bus or train in our 

daily commute, and often walk to connect between networks or to our final destination. 

A transportation network consists of a set of n non-intersecting roads, where each road 
has a speed. Thus a transportation network is usually modelled as a plane graph T{S,C) in 
ff^ the Euclidean plane (or some other metric) whose vertices S are nodes and whose edges C are 

roads. Furthermore, each edge has a weight a S (0, 1] assigned to it. One can access or leave a 
road through any point on the road. In the presence of a transportation network, the distance 
^ between two points is defined to be the shortest elapsed time among all possible paths joining 

the two points using the roads of the network. The induced distance, called d-j-, is called a 
transportation distance. 
^> Using these notations the problem at hand is as follows: 

X 





Fig. 1. Illustrating a quickest path from a source point s to a destination point t. 
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Problem 1. Given two points s and t in and a transportation network T{S,C) in the 
Euclidean plane. The problem is to find a path with the smallest transportation distance 
from s to t, as shown in Fig. 1. 

Most of the previous research has focussed on shortest paths and Voroni diagrams. Abel- 
lanas et al. [1,2] started work in this area considering the Voronoi diagram of a point set 
and shortest paths given a horizontal highway under the Li-metric and the Euclidean metric. 
Aichholzer et al. [4] introduced the city metric induced by the Li-metric and a highway net- 
work that consists of a number of axis-parallel line segments. They gave an efficient algorithm 
for constructing the Voronoi diagram and a quickest path map for a set of points given the 
city metric. Gorke et al. [14] and Bae et al. [5] improved and generalised these results. 

In the case when the edges can have arbitrary orientation and speed, Bae et al. [6] presented 
algorithms that compute the Voronoi diagram and shortest paths. They gave an algorithm 
for Problem 1 that uses 0{n^) time and O(n^) space. This result was recently extended to 
more general metrics including asymmetric convex distance functions [7] . 

In this paper we improve on the results by Bae et al. [6] and give an 0{n^ log n) time 
algorithm using O(n^) space. Furthermore, we introduce the (approximate) query version. 
That is, a transportation network with n roads in the Euclidean plane can be preprocessed 
in 0(n^ log n) time into a data structure of size 0{n^/s'^) such that given any two points s 
and t in the plane a (1 + £)-approximate quickest path between s and t can be answered in 
0(l/e^ • logn) time. For the query structure we assume that the minimum and maximum 
weights, ctmin and ctmax, are constants in the interval (0, 1] independent of n, the exact bound 
is stated in Theorem 8. 

This paper is organized as follows. Next we prove three fundamental properties about an 
optimal path among a set of roads. Then, in Section 3, we show how we can use these properties 
to build a graph of size O(n^) that models the transportation network, a source point s, a 
destination point t and, contains a quickest path between s and t. In Section 4 we consider 
the query version of the problem. That is, preprocess the input such that an approximate 
quickest path query between two query points s and t can be answered efficiently. Finally, we 
conclude with some remarks and open problems. 

2 Three Properties of an Optimal Path 

In this section we are going to prove three important properties of an optimal path. The 
properties will be used repeatedly in the construction of the algorithms in Sections 3 and 4. 
Consider an optimal path V and let Cp = (ci, C2, . . . , c^) be the sequence of roads that are 
encountered and followed in order along V, that is, the sequence of roads on which the path 
changes direction. For example, in Fig. 1 that sequence would be {ci,C4,cq) but not include 
C3 since the path does not follow C3. For any path Vi, let wt(Vi) denote the cost of the path 
Vi. Note that a road can be encountered, and followed, several times by a path. For each 
road Ci G C-p let Si and ti be the first and last point on Cj encountered for each occasion by V. 
Without loss of generality, we assume that tj+i lies below Sj when studying two consecutive 
roads in C-p. We have: 

Property 1. For any two consecutive roads Cj and c^+i in Cp the subpath of V between ti and 
Sj+i is the straight line segment 



The endpoints of a road Cj = {uj,Vj) are denoted start point (uj) and end point (vj), as 
they appear along the road direction. 

Consider a segment connecting two consecutive roads Cj and Cj+i along C-p. Let 

<pi+i denote the angle Z{ti, Si+i,Ui+i), as illustrated in Fig. 2(a). 

Property 2. If s^+i lies in the interior of Cj+i then ^j+i = arccos(ai+i). 

Proof. For simplicity rotate C such that Cj+i is horizontal and lies below tj, as shown in 
Fig. 2(a). Let r denote the orthogonal projection of ti onto the ray containing Cj+i (not 
necessarily on Cj+i and let h = \tir\. We have: 

i, I h Jill. cos(/)i+i 

ijSi+i = -. — : and rs^+i = h ■ — — . 

sm (pi+i sm <pi+i 

Thus, the cost of the path from ti to ij+i along C-p as a function of 0^+1 is: 

JWi+l) = \tiSi+l \ + Oj+l • Si+ltj+l = h Oi+l • ( rtj+i - /l • ). 

sm^j+1 sm(^i+i 
Differentiating the above function with respect to (f)i+i gives: 

/i(cOS(/)i+l - Oj+l) 



COs2 - 1 



Setting f'{(pi^i) = the resulting function gives that the minimum weight path between ti 
and fj+i along C{V) is obtained when 

(f)i+i = arccos(ai+i). 



□ 



Property 3. There exists an optimal path V' of total cost wt{V) with C-p' = C-p that fulfills 
Properties 1-2 such that for any two consecutive roads Cj and q+i in Cp' the straight-line 
segment {ti, Sj+i) of V' must have an endpoint at an endpoint of Ci or Cj+i, respectively. 

Proof. Assume the opposite, i.e., {ti, s^+i) does not coincide with any of the endpoints of Cj or 
Ci+i. Consider the three segment path from Sj to tj+i, that is, (sj, tj), {ti, Sj+i) and (sj+i, ij+i). 
The length of this path is: 

C^i ■ I •Si, ill + l^i, Si+ll + Cki+l • 

According to Lemma 2 the orientation of (tj,Sj+i) is fixed, which implies that the weight of 
the path is a linear function only depending on the position of ti (or Sj+i). Hence, moving ti in 
one direction will monotonically increase the weight of the path until one of two cases occur: 
(1) either ti or Sj+i encounters an endpoint of Cj or Cj+i, or (2) ti = Si or Sj+i = t^+i. If (1) 
then we have a contradiction since we assumed {ti, Sj+i) did not coincide with any endpoint. 
And if (2) then we have a contradiction since V must follow both Cj and Cj+i (again from the 
definition of Cp). □ 




Fig. 2. (a) Illustrating Property 2. (b) Defining PRf{s,Ci) and (c) PRb{t,Ci). 



3 The basic case 

In this section wc consider Problem 1, that is, as input we are given a source point s, a 
destination point t and a transportation network T{S,C) and the aim is to find a path with 
minimum transportation distance from s to t. Our algorithm will construct a graph G that 
models the set C of roads and quickest paths between the roads. Using the three properties 
shown in the previous section we will show that if a shortest path in G between s and t has 
cost W then an optimal path between s and t has cost W. The optimal path can then be 
found by running Dijkstra's algorithm [13] on G. 

Fix an optimal path V that fulfills Properties 1-3 (we know such a path exists). If V 
follows Ci and Cj+i then the path between Cj and Cj_|-i (a) is a straight line segment, (b) the 
segment (ti, Sj+i) forms a fixed angle with Cj_|-i, and (c) at least one of its endpoints coincides 
with an endpoint of q or Cj+i. These three properties suggest that V has a very restricted 
structure which we will try to take advantage of. 

Let PRf(p, Ci) be the projection of a point p onto a road Cj, if it exists, such that the angle 
Z(|), PRf{p, Ci),Ui) is as shown in Fig. 2(b). Furthermore, let PRh{p, a) be the projection 
of point p on a road Cj, if it exists, such that the angle Z(j), PRi,{p, Ci),Vi) is (pi, as shown in 
Fig. 2(c). 

Consider the graph G{V, E) with vertex set V and edge set E. Initially V and E are 
empty. The graph G is defined as follows: 

1. Add s and t as vertices to V . 

2. Add the nodes in S as vertices to V . 

3. For every road Cj G C add the point PRf(s, q) (if it exists) as a vertex to V and add the 
directed edge {s, PRf{s,Ci)) (if it exists) with weight \sPRf{s,Ci)\ to E, see Fig. 3(a). 

4. For every road Ci £ C add the point PRi,{t, Ci) (if it exists) as a vertex to V and add the 
directed edge {PRb{s,Ci),t) (if it exists) with weight \PRi,(s, Ci)t\ to E. 

5. For every pair of roads Ci,Cj G C add the following points (if they exist) as vertices to 
V: PRf{vi,Cj), PRf{ui,Cj), PRf{uj,Ci) and PRf{vj,Ci). Add the following four directed 
edges to E (if their endpoints exist): {vi, PRf{vi,Cj)), {ui, PRf{ui, Cj)), {vj, PRf{vj,Ci)) 
and {uj, PRf{uj,Ci)). The weight of an edge is equal to the Euclidean distance between 
its endpoints, see Fig. 3(b). 

6. For every pair of roads Ci,Cj G C add the directed edges {vi,Uj) with weight \viUj\ and 
{vj,Ui) with weight \vjUi\ to E. 

7. For every road q consider the vertices of V that correspond to points on in order from 
Ui to Vi- For every consecutive pair of vertices Xj,Xj+i along q add a directed edge from 
Xj to Xj+i of weight CKj • as shown in Fig. 3(c). 



Fig. 3. Illustrating how the graph is built. 



Lemma 1. The graph G contains O(n^) vertices and O(n^) edges and can be constructed in 
time 0{n^ logn). 

Proof. For every pair of roads we construct a constant number of vertices and edges that 
are added to V and E, thus O(n^) vertices and edges in total. For the first five steps of the 
construction the time to construct the vertices and edges is hnear with respect to the size of 
the graph, since every edge and vertex can be constructed in constant time. In step 6 we need 
to sort 0{n) vertices along each road, thus 0(n^ logn) time in total. □ 

The following observation follows immediately from the construction of the graph and 
Properties 1-3. 

Observation 1 The shortest path between s and t inG has cost W if and only if the minimum 
transportation distance from s to t has cost W . 

By simply running Dijkstra's algorithm [13], implemented using Fibonacci heaps, on G 
gives the main result of this section. 

Theorem 1. A path with minimum transportation distance between s and t can be computed 
in O(n^logn) time using 0{n^) space. 

3.1 Shortest paths among polygon obstacles 

In this section we briefly discuss how the above algorithm can be generalised to the case 
when the plane contains polygonal obstacles. As input we are given a source point s, a 
destination point t, a transportation network T{S,C) in the Euclidean plane and a set O of 
k non-intersecting obstacles of total complexity m. 

Every edge of an obstacle can be viewed as an undirected road (or two directed edges) 
with associated cost function a = 1. Consequently, the edge connecting a road and an obstacle 
edge along the optimal path has the three properties described in Section 2. However, while 
constructing the graph we have to add one additional constraint, namely, no edge in E can 
intersect an obstacle. According to these properties we are going to build the graph G that 
models the set of roads and obstacles. 

There are several methods to check if a segment intersects an obstacle [3, 12, 15]. We will 
use the data structure by Agarwal and Sharir [3] which has 0[m^^^ j \flJ) query time using 
©(L^"*"^) preprocessing and space for m < L < rr? . Using this structure with L = vr? gives 
us the following results: 



Lemma 2. The graph G has size 0{N'^~^^) and can be constructed in time 0{N'^~^^), where 
N = n + m. 

By simply running Dijkstra's algorithm, implemented using Fibonacci heaps, on G gives: 

Theorem 2. A collision-free path with minimum transportation distance between s and t 
among O can be computed in 0(N'^ logN) time using 0{N'^) space, where N = n + m. 

4 Shortest Path Queries 

In this section we turn our attention to the query version. That is, preprocess C such that 
given any two points s and t in find a quickest path between s and t among C effectively. 
We will present a data structure P that returns an approximate quickest path. That is, given 
two query points s and t, and a positive real value e, the data structure returns a path 
between s and t having transportation distance at most (1 + e) times the cost of an optimal 
path between s and t. 

To simplify the description we will start (Sections 4.1-4.2) with the case when t is already 
known in advance and we are only given the start point s as a query. Then in Section 4.4 we 
generalize this to the case when both s and t are given as a query, and in Section 4.5 we show 
how one can improve the preprocessing time and space requirement using the well-separated 
pair decomposition. 

Fix an optimal path V that fulfills Properties 1-3 (we know such a path exists), and 
consider the first segment I = (s,si) of V. Obviously si must be either a start/endpoint of 
a road (type 1) or an interior point of a road q G C (type 2) such that i and q form an 
angle of (f)i (ignoring the trivial case when si = t). We will use this observation to develop an 
approximation algorithm. The idea is simple. Build a graph G{V, E), as described in Section 3, 
with T and t as input. Compute the cost of the quickest path between t and every vertex in 
V . Now, given a query point s, find a good vertex si in V to connect s (either directly or via 
a 2-link path) to and then lookup the cost of the quickest path from si to t in G. Obviously 
the problem is how to find a "good" vertex. In the next subsection we will select a constant 
number of candidate vertices such that we can guarantee that at least one of the vertices will 
be a "good" candidate, i.e., there exists a path, that fulfills Properties 1-3, from s to t via 
this vertex that has cost at most (1 + e) times the cost of an optimal path. 

4.1 Finding good candidates: Type 1 and Type 2 

Let Sc denote the set of the endpoints (both start and end) of the roads in C and let s be the 
query point. As described above we will have two types to consider, and we will construct a 

set Vi for the type 1 cases and a set 2^2 for the type 2 cases. The first set, Pi, is a subset of Sq 
and the second set, T)2, '\s a, set of 3-tuples that will be used by the query process (described 
in Section 4.2) to calculate the quickest path. 

Type 1: For the point set T>i we will use the same idea as is used in the construction of 
^-graphs [17]. Partition the plane into a set oi k = max{9, ^} equal size cones, denoted 
Xi, . . . with apex at s and spanning an angle of 6 = 2TT/k, as shown in Fig. 4a. For 
each cone X the set Pi contains a point r, where r is a point in 5c H X whose orthogonal 
projection onto the bisector of X is closest to s. The following holds: 



Lemma 3. Given a point set Sc and a positive constant e one can preprocess Sc into a data 
structure of size 0{n/e) in 0{l/e • nlogn) time such that given a query point s the point set 
Vi, of size at most 36Tr/e, can be reported in 0{l/s ■ logn) time. 

Proof. Given a direction d and a point s let £d{s) denote the infinite ray originating at s 
with direction d, see Fig. 4b. Let C{s,d,9) be the cone with apex at s, bisector £d{s) and 
angle 9. It has been shown (see for example Section 4.1.2 in [17] or Lemma 2 in [10]) that Sc 
can be preprocessed in O (nlogn) time into a data structure of size 0(n) such that given a 
query point s in the plane the data structure returns the point in C(s, d, 9) whose orthogonal 
projection onto ld{s) is closest to s in O(logn) time. We have 367r/£ directions, thus the 
lemma follows. □ 



Fig. 4. (a) Partitioning the plane into k cones, (b) Selecting the point whose orthogonal projection onto the 
bisector of X. 



Type 2: It remains to construct the set T>2 of 3-tuples. Unfortunately the construction 
might look unnecessarily complicated but hopefully it will become clear, when we prove the 
approximation bound (Section 4.3) why we need this construction. Before constructing T>2 we 
need some basic definitions. 

Let 

Q^max and ctmin be the maximum and minimum weight of the roads in C. Partition the 
set C of roads into a minimum number of sets, Ci, . . . , Cm, such that the orientation of a road 
in Ci is in [(i — 1)^ • amin, i9 ■ CKmin), where 9 = e/18. Partition each set Cj, 1 < i < m, into 6 
sets Ci^i, . . . ,Ci,b such that the weight of a road in C^j is in [amin • (1 + e)-'"^, amin • (1 + £)■')• 

For every i, j, 1 < i < m and 1 < j < b, define two directions 7up(i,i) and 'Jdownihj) as 
follows (see also Fig. 5a). Consider an infinite ray T with orientation (i — 1)9 having weight 
" = «min • For simplicity we rotate the ray such that it is horizontal directed from left 
to right. Let jupihj) be the direction of a ray r originating from below T such that r and T 
meet at an angle of arccos(a). The direction "Jdownihj) is defined symmetrically but with a 
ray originating from above T. 

Given a point p on a road c let N(p, c) denote the nearest vertex of G to p along c. Note 
that N{p, c) must lie between p and the end point of c. 

Now, we are ready to construct 2^2 • When given the query point s construct the set P2 as 
follows. For each i,j, 1 < i < m and 1 < j < b, shoot a ray r^p originating from s in direction 
lup{i-,j) and one ray rdown in direction Jdownih j)- If i"up hits a road in Cij then let c„p be 
the first road hit and let Pup be the point hit on c„p, as illustrated in Fig. 5b. The 3-tuple 




[cup,Pup, ^iPup, Cup)] is added to V2. If r^own hits a road in Cjj then let Cdown be the first 
road hit and let pdown be the point hit on Cdown- The 3-tuple [cdown,Pdown, N{pdown, Cdown)] is 
added to V2. 

Lemma 4. Given a transportation network T{S,C) with n roads in the Euclidean plane and 
a positive constant s one can preprocess T in O(nlogn) time into a data structure of size 
0(n) such that given a query point s the set D2 can be reported in 0( — - — 2 dognlogi , ^ ^max) 

time. The number of 3-tuples in P2 is 0( — ^— ^ • logi ^™'). 

Proof. The preprocessing consists of two steps: (1) partitioning C into the sets Cjj, 1 < i < m 
and 1 < j < b, and (2) preprocessing each set Cij into a data structure that answers ray 
shooting queries efficiently. 

The first part is easily done in 0(ri log n) time by sorting the roads first with respect to 
their orientation and then with respect to their weight. 

The second step of the preprocessing can be done by building two trapezoidal maps 
Tup{Cij) and Tdown{Ci,j) (also know as a vertical decomposition) of each set Cij as follows 
(see Fig. 5c). Rotate Cij such that ^upihj) is vertical and upward. Build a trapezoidal map 
Tup{Ci,j) of dj as described in Chapter 6.1 in [8]. Then preprocess Tup{Cij) to allow for planar 
point location. Note that every face in T{C) either is a triangle or a trapezoid, and the left 
and right edges of each face (if they exists) are vertical. The trapezoidal map Tdowni^ij) can 
be computed in the same way by rotating Cij such that 'Ydownihj) is vertical and upward. 
The total time needed for this step is O(nlogn) and it requires 0(n) space. 

When a query point s is given, two ray shooting queries are performed for each set Cij. 
However, instead we perform a point location in the trapezoidal maps Tup{Cij) and Tdown{Cij). 
Consider Tup{Cij) and let / be the face in the map containing s. The top edge of / corresponds 
to the first road Cup hit by a ray emanating from s in direction 7Mp(i,i). When Cup is found we 
just add to 2^2 the first vertex on Cup in G to the right of s. The same process is repeated for 
Tdown{Cij)- Performing the point location requires O(logn) time per trapezoidal map, thus 
the total query time is 0(m6 logn) = 0{ .\ 2 • lognlogi , ^.(amax/amin))- □ 




Fig. 5. (a) Illustrating the definition of 'yupihj) and 'ydo-wn{i, j)- (b) Pup is the first point hit by the ray and 
N{cup,Pup) is the neaxest neighbour oi pup along Cup- (c) The trapezoidal map of a set Cij and the query 
point s. 



4.2 The preprocessing and the query 



In this section we will present the remaining data structure, define the preprocessing and 
show how a query is processed. 

Preprocessing In Section 3 we showed how to build a graph given two points and the 
transportation network T{S, C). The first step of the preprocessing is to build a graph G(V, E) 

of C and the destination point t (without the source point). Next compute the shortest path 
from every vertex in G to the vertex corresponding to t. Since the complexity of G is quadratic, 
this step can be done in 0(n^ log n) time and 0{'n?') space using Dijkstra's algorithm^. The 
distances are saved in a matrix M. Finally, we combine the above results with Lemmas 3 
and 4 and get: 

Theorem 3. The preprocessing requires 0{v?\ogn) time and 0{n'^) space. 

Query As a query we are given a point s in the plane. First we compute the two sets Pi 

and For each vertex si in T>i compute the quickest path via si, that is, \ssi\ + M\si,t]. 
The quickest among these paths is denoted Vi- The 3-tuples in T>2 require slightly more 
computation. For each 3-tuple [c, si, S2] consider the path from s to t using (s, s\) and (si, S2). 
The cost of the path can be calculated as |ssi| + a ■ \s\S2\ + M[s2,i], where a is the weight 
of road c. The quickest among these paths is denoted V2- 

The quickest path among Vi, V2 and the direct path from s to i is then reported. 

Theorem 4. A query can he answered in time ^ ^2 • lognlog^^^ ^^)- 

Proof. As above we divide the analysis into two parts: type 1 and type 2. 

Type 1: According to Lemma 3 the number of type 1 candidate points is at most 2>1-k je and 
can be computed in time 0(l/e • logn). Computing the cost of a quickest path for a point 
in Pi can be done in constant time. Thus, the query time is 0(l/£:logn). 

Type 2: According to Lemma 4 the number of 3-tuples in P2 is 0( — ^-^^ • logi , , amax/«min) 
and can be computed in time ^ ,2 ' log7ilogi_|_£ amax/ctmin)- Each clement in P2 is 

then processed in time, thus 0( — ^—2 • logralogi,, amax/ctmin) in total. 

^min^ 

Summing up we get the bound stated in the theorem. □ 
4.3 Approximation bound 

Consider an optimal path V and let i = {s, si) be the first segment of V. According to 
Property 1 it is a straight-line segment. For any two points pi and p2 along V, let 6GiPijP2) 
denote the cost of V from pi to p2. Define s(i) to be the sector with apex at s, bisector on i, 
interior angle k = e/18 and radius (1 -|- Ak) ■ \ssi\, as shown in Fig. 6(a). 

Lemma 5. Let e be a positive constant. If s{£) CiSc 7^ then: 

wt{V) < wtiVi) <{! + £)■ wt{V). 
^ The SSSP has to be performed from i on a graph G' where every edge has swapped direction. 



Proof. The proof will be shown in two steps: (1) first we prove that for every cndpoint p 
within s{i) the quickest path, denoted V{p), from s to t using the segment {s,p) has cost at 
most {1 + 12k) ■wt{V), and then (2) we prove that wt{Vi) has cost at most {1 + 3k) ■wt{V{p)). 
Combining the two parts proves the lemma. 

Part 1: Consider any point p within s{£), and let 'P{p) denote the quickest path from s to 
t using the segment {s,p). We have: 



Part 2: As above let p be any endpoint of C within s{i), and let V{p) be the quickest path 
from s to t using the segment \sp\. 

Consider the set X>i as described in Section 4.1 and assume without loss of generality that 
p lies in a cone X. If there exists a point q in T>i such that q = p then we are done since 



Otherwise, there must exist another point q in Vi such that q E X and whose orthogonal 
projection q' onto the bisector of X is closer to s than the orthogonal projection p' of p onto 
the bisector of X. 



wt{V{p)) = \sp\ + 6G{p,t) 
< \sp\ + \psi \ 



< (1 + Ak) ■ \ssi \ + Y (cos - - 

< (1 + 12k) ■ \ssi\ + 6G{si,t) 

< {l + 12K)-SG{s,t) 




ssi\ + SG{si,t) 



sq\+SG{q,t)=wt{V{p)). 



wt{V{q)) < \sq\ + \qp\ + 5G(p,t) 



< \sq\ + \qq'\ + \q'p'\ + \p'p\ + SG{p,t) 



< \sQ\i 1 ) + kVl + \sp\ ■ smK + 6Gip,t) 



cos K cos K 



< {\sq'\ + \q'p'\){ \ ) + \sp\ ■ smK + SG{p,t) 



cos K COS K 




< {1 + 2 sin K) + 6G{p,t) 



cos K 



< |sp| • (1 + 3cotK) + 5G{p,t) 

< (1 + 3k) ■ wt{P{p)) 



In the last step we used that k = e/18 < 2-!t/9. 
Now we can combine the two results as follows. 



wt{Vi) < wt{V{q)) 



< {l + 12K)-SG{p,t) 

< {1 + 3k) ■ {l + 12K)-6G{s,t) 
<{l + s)-6Gis,t) 



(from Part 1) 
(from Part 2) 

(since k = e/18 and £ < 1) 



This completes the proof of the lemma. 



□ 




Fig. 6. (a) Illustrating the definition of s{£). (b) Illustrating the setting in Lemma 6. 

Lemma 6. Let s < 1 be positive constants. If s(C) n<Sc = then: 

wt{r) < wt{V2) <(! + £)• wt{P). 

Proof. As above let (s, si) be the first segment of V, where si hes on a road ci. Assume w.l.o.g. 

that ci is belongs to the set of roads Cij as defined in Section 4.1 (Type 2). Rotate the input 
such that ci is horizontal, below s and going from right to left. Consider the construction of 
the 3-tuples in Type 2, and let [c'i,s'^,p] be the 3-tuple reported when Cij was processed in 
the direction towards ci. See Fig. 7(b). 

Consider the shortest path from s to t using the segment {s,s'i). We will have two cases 
depending on s[: either (1) si and s'^ lie on the same road, or (2) they lie on different roads. 

Case 1: If s[ and si both lie on road ci (with cost function ai) then we have: 

wt{V2) < \ss[\+ai • \s[,ti\ +5G{ti,t) 

< \ss'i\ + ai ■ \s[si\ + ai ■ \siti\ + dG{ti,t) 

< \ss'i\ + ai ■ \s[si\ + 5Gisi,t) 

< (1 + 4k) ■ \ssi \ + 5k ■ \ssi \ + daisi, t) (since 6 < k) 

< (1 + 9k) • |ssi| +dG{si,t) 

<{l + <dK)-5G{s,t) 

<{l + e)-6G{s,t) 

That completes the first part. 

Part 2: If s'^ and si lie on different roads c'^ and ci, respectively, then c'^ must lie between 
s and ci. This follows from the fact that s{tj does not contain any endpoints and c/^ is the 
first road hit. Furthermore, there must exist a an edge (g', q) e E such that q' lies on to 
the left of s[ and q lies on ci to the right of ti and Z{ui,q,q') = arccosai. See Fig. 7(b) for 
an illustration of case 2(a). 

Consider the situation as depicted in Fig. 7(b). We will prove that the cost of the path 
from s to ti via s'^ , q' and q is almost the same as the cost of the optimal path from s to ti 



(that goes via si). Recall that the cost function of ci and is cti and the cost function of c'^ is 
a[. Furthermore, let r be the intersection point between cf^ and (s, si). 

Note that the cost of the path from r to g via q' is maximized if rq' forms an angle of with 
the horizontal line and q' lies above r, as shown in Fig. 7(b). Furthermore, \qq'\ = \q'p\ + \pq\ 
and \pq\ = \rsi\. 

a'l ■ \rq'\ + \q'p\ < a\ • \pr\ • cos 9 + \pr\ ■ sin9 

< \pr\ia'i + s- amin/18) 

< a'l ■ \pr\{l + e) 

Putting together the bounds we get: 

wt{V2) < \ss'i \ +ai ■ \ s[q'\ + \q'q\ + 5G{q,t) 

< \ss'i\ + a'l ■ \s'ir\ + a'l ■ \rq'\ + {\q'p\ + \pq\) + 5G{q,t) 

< (1 + 9k) • \sr\ + a'l • \rq'\ + \q'p\ + |rsi| + dciq, t) (see Part 1) 

< (1 + 9k) ■ \sr\ + |rsi| + a'l ■ \pr\ (1 + e) + 6Giq, t) 

< (1 + e) • {\sr\ + \rsi\ + a'l ■ \siq\) + fc(g, t) 

< {l + e)-wt{P) 

This completes the proof of Lemma 6. □ 




Fig. 7. Illustrating (a) Case 1 and Case 2 (b) in the proof of Lemma 6. 



We can summarize this section (Lemmas 5 and 6, Theorems 3 and 4) with the following 
theorem: 



Theorem 5. Given a transportation network T with n roads in the Euclidean plane, a des- 
tination point t and a positive constant e, one can preprocess T and t in O(n^logn) time 



and 0{n?) space such that given a query point s, a {1 + e)- approximate quickest path can be 
calculated in 0( — ^ • lognlogi , , "wx-j ^j^g, 

4.4 General case 

In this section we turn our attention to the query version when we are given two query 

points ,s and t in and our goal is to find a quickest path between s and t among C. The 
idea is tlie same as in the previous section. That is we perform the exact same preprocessing 
steps as in the previous section (omitting the destination point t) , but with the exception that 
M contains all-pair shortest distances. Using Johnson's algorithm [16] the all-pairs shortest 
paths can be computed in 0(n^ log n) using O(ra^) space. 

Given a query we compute the two candidate sets of type 1 and type 2 for both s and t. 
For each pair of elements p G T^i{s) U 'D2{s) and q £ T^i{t) U T>2{t) compute a path between s 
and t as follows: 

if p GVi and q GVi then \sp\ + M\p, q] + \qt\ 

if p G Pi and q = [c, si, S2] G ^2 then \sp\ + M\p, S2] + a(c) • |s2Si| -|- |sis| 
if p = [c, si, S2] G ^2 and q & Vi then |ssi| + a(c) • |siS2| + M[s2,q] + \qt\ 
iip= [ci, 81,82] G 1)2 and q = [c2,ti,t2\ G T>2 then |ssi| + a{ci) ■ \siS2\ + M[s2,t2] + a{c2) ■ 

\t2tl\ + \tlt\ 

Note that for a specific pair p, q the calculated distance might not be a good approximation 
of the actual distance. However, for the shortest path there exists a pair p, q such that the 
distance is a good approximation. 

By putting together the results, we obtain the following theorem: 

Theorem 6. Given a transportation network T with n roads in the Euclidean plane and a 
positive constant e, one can preprocess T in 0(n^ log n) time using O(n^) space such that 
given two query points s and t a (l + e)- approximate quickest path between s and t can be 
calculated in 0(( — - — o • logi4_=- ^^^Y ■ logn) time. 

4.5 Improving the complexity using the well-separated pair decomposition 

The bottleneck of the preprocessing algorithm is the fact that shortest paths are computed 
in a graph of quadratic complexity. Is there a way to get around this? Since it suffices to 
approximate the shortest paths we can reduce the number of shortest path queries from O(n^) 
to O(n^), that is, linear in the number of vertices, using the well-separated pair-decomposition 
(WSPD). 

Definition 1 ([11]). Let s > be a real number, and let A and B be two finite sets of points 

in W^. We say tha,t A and B are well-separated with respect to s, if there are two disjoint d- 
dimensional halls Ca and Cb, having the same radius, such that (i) Ca contains the hounding 
box R{A) of A, (i) Cb contains the bounding box R{B) of B, and (ii) the minimum distance 
between Ca and Cb is at least s times the radius of Ca- 

The parameter s will be referred to as the separation constant. The next lemma follows 
easily from Definition 1. 



Lemma 7 ([11]). Let A and B be two finite sets of points that are well- separated w.r.t. s, 
let X and p be points of A, and let y and q be points of B. Then (i) \xy\ < (1 + 4/s) • \pq\, and 
(ii)\px\<{2/s)-\pq\. 

Definition 2 ([11]). Let S be a set of n points in W^, and let s > be a real number. A 
well-separated pair decomposition (WSPD) for S with respect to s is a sequence of pairs of 
non-empty subsets of S, {Ai,Bi), . . . , {Am,Bm), such that 

1. Ai{~\Bi = 0, for all i = 1, . . . ,m, 

2. for any two distinct points p and q of S, there is exactly one pair {Ai, Bi) in the sequence, 

such that (i) p G Ai and q £ Bi, or (ii) q £ Ai and p G Bi, 

3. Ai and Bi are well-separated w.r.t. s, for 1 < i < m. 

The integer m is called the size of the WSPD. 

Callahan and Kosaraju showed that a WSPD of size m = 0{s^n) can be computed in 
0{s'^n + nlogn) time. 

Construct the graph G{V, E) of T, as defined in Section 3. Compute a WSPD {{Ai, Bi)}^^-^ 
of V with respect to a separation constant s = — - — , where r < 1 is a positive constant 
given as input. Then for each well-separated pair {Ai,Bi) pick two arbitrary points a £ Ai 
and b £ Bi as representative points, and calculate the shortest path in G between a and b. 
All paths are stored in a matrix M'. According to Definition 2, we have O(n^) well separated 
pairs. It follows that the number of the shortest path queries in M is O(n^). 

The queries are performed in almost the same way as above. The only difference is how 
the cost of the path between two points is calculated. Assume that we want the minimum 
transportation distance between two query points p and q. According to Definition 2 there 
exists a well-separated pair {A, B} such that p £ A and q £ B, or vice versa. Furthermore, let 

and r^ be the representative points of A and B, respectively. Instead of using the value 
in M\p, q] we approximate the cost by \prA\ + M'[rA, tb] + \rBQ\- 

Theorem 7. 

^g{p, q) < \prA\ + M'[rA, rs] + {rsql < (1 + r) • Sg{p, q). 
Proof. The left inequality is immediate, thus we will focus on the right inequality. 

\prA\ +M'[rA,rB] + \rBq\ = \prA\ +<^g(m,?^b) + \rBq\ 

< 2\prA\ + 2\rBq\ + Scip, q) 
<4-2/s-\pq\+6Gip,q) 

< r ■ amin • \pq\ + SG{p,q) 

< il + T)-5G{p,q) 

Where the last step follows by using \pq\ < Sg{p, q) / amin. D 

Given a positive constant cr < 1 we can obtain the claimed bounds by setting the constants 
appropriately (for example r = a/3 and e = a/3). We get: 

Theorem 8. Given a transportation network T{S, C) with n roads and a positive constant e, 

one can preprocess C in 0(( — ^^)'^logn) time using 0(( — - — )^) space such that given two 
query points s and t a (1 + e)- approximate quickest path between s and t can be calculated in 
0{{^^ • logi+^ 1^)2 • logn) time. 

Assuming amin and ctmax being constants the bounds can be rewritten as: preprocessing 
is 0(n^ logn), space is 0{n^ /e^) and the query time is 0(l/£:^ • logn). 



5 Concluding Remarks 



We considered the problem of computing a quickest path in a transportation network. In the 
static case our algorithm has a running time of 0(n^ log n) which is a linear factor better 
than the best algorithm known so far. We also introduce the query version of the problem. 

There are many open problems remaining. For example, can one develop a more efficient 
data structure that has a smaller dependency on amin and e? Also, is there a subquadratic 
time algorithm for the static case? If not, can we prove that the problem is 3sum-hard? Bae 
and Chwa [6] proved that the shortest path map for a transportation metric can have f2('n?) 
size, which indicates that the problem may indeed be 3sum-hard. 
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